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We study the depinning transition for models representative of each of the two universahty classes 
of interface roughening with quenched disorder. For one of the universality classes, the roughness 
exponent changes value at the transition, while the dynamical exponent remains unchanged. We 
also find that the prefactor of the width scales with the driving force. We propose several scaling 
relations connecting the values of the exponents on both sides of the transition, and discuss some 
experimental results in light of these findings. 

PACS numbers: 47.55.Mh 68.35.Fx 
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Interface roughening in the presence of quenched dis- 
order has recently been the focus of great interest [Q. 
In the typical d-dimensional interface described 

by a height h{x.,t) moves in a (d-l- l)-dimensional disor- 
dered medium. The randomness of the medium can be 
described by a quenched disorder ri{x, h). 

The continual motion of the interface in the medium 
requires the application of a driving force F . For driving 
forces smaller than a critical value Fc, the interface be- 
comes pinned by the disorder after some finite time. For 
F > Fc, the interface moves with a constant velocity v. 
Hence, the velocity can be identified as the order param- 
eter of the depinning transition. Close to the transition, 
the velocity of the interface scales as ^ where 9 
is the velocity exponent, and f = {F — Fc)/Fc is the 
reduced force. 

Experimental studies of the scaling of the local inter- 
face width, w , where t is the window of observation, 
reveal roughness exponents a in the range 0.6 — 1.0 [^,^. 
Several models in which quenched disorder plays an es- 
sential role have been proposed with the goal of explain- 
ing the experimental results (for reviews, see e.g. jl]). For 
these models, both the scaling of w with I and the scaling 
of the global interface width W with the system size L, 
were studied. 

A numerical study Q], recently confirmed by analyti- 
cal arguments revealed that most of the models in- 
troduced so far could be organized into two universal- 
ity classes. A group of these models can be mapped, at 
the depinning transition, onto directed percolation (DP) 
for d = 1, or directed surfaces (DS) for d > 1 [^,0. 
These models are referred to as directed percolation de- 
pinning (DPD) models. The mapping allows the deter- 
mination of a from the correlation length exponents of 
DP or DS. Reference |Q showed that the DPD models 
can be described by a stochastic differential equation of 
the Kardar-Parisi-Zhang (KPZ) type with quenched 
disorder 101 



dh 
'dt 



F + V^h + X(yhf +r^{x, h), 



(1) 



the depinning transition. 

A number of different models, belonging to the sec- 
ond universality class, to which we refer to as quenched 
Edwards-Wilkinson (QEW), were found to have either 
A = or A ^ at the depinning transition Q. At the 
depinning, they can be described by the EW equation 
]lO| with quenched disorder h% 



dh 
'dt 



F + V2/i + ry(x, h). 



(2) 



This equation has been studied by means of the func- 
tional renormalization group |12|. Numerical studies of 
the pinned phase yield a ~ 0.97 ]13|] and a ~ 1.25 ||l^ in 
(1-1-1) dimensions. 

A remaining problem is the connection of these univer- 
sality classes to the experiments. We observe that most 
numerical and analytical studies focus on the "pinned 
phase" {F < Fc) while nearly all experimental results 
are for the "moving phase" {F > Fc). 

In this Letter, we study the moving phase in (1 + 1)- 
dimensions for models representative of each of the two 
universality classes. For the DPD universality class, we 
find that both a and the growth exponent (3 — defined 
hy W ^ t^ for t ^tx {L), where tx {L) is the saturation 
time — have different values on both sides of the de- 
pinning transition. However, we find that the dynamical 
exponent z — a/ (3 has the same value in both phases. 

First, we consider the QEW universality class. In 
(1 + l)-dimensions, we study a Hamiltonian model de- 
fined as 



^ ^ H [(/ifc+i - /ifc)' - Fhk + ri{k, hk 



(3) 



fc=i 



where the coefficient A of the nonlinear term diverges at 



Here, the first term represents the elastic energy that 
tends to smooth the interface, and r]{k, hk) is an uncorre- 
lated random number, which mimics a random potential 
due to the disorder of the medium. In the simulation, a 
column k is chosen and its height is updated to {hk + 1) if 
the change in (|^) is negative. Thus, only motions that de- 
crease the total energy of the system are accepted. Back- 
ward motions are neglected since these are rare events. 
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We start by studying the scaling of the local width w 
in a window of size £, for different driving forces. In Fig. 
|l|(a), we show the local width for the pinned and the mov- 
ing phases. Consecutive slopes yield roughness exponents 
Q!p ~ 0.92 for the pinned phase, and am — 0.92 for the 
moving phase (Table ^. We also detect a dependence of 
the local width on L, that can be described by a scaling 
of the form ifisi 



For the pinned phase we have, according to 



(4) 



where is a scaling function that for w <C 1 scales as 
u"". We find ~ 1.23. 

We study the scaling of W with t, and find f3p ~ 0.85 
and Prri — 0.86. These results imply that z remains un- 
changed at the transition. 

The study of the scaling of the local width for the mov- 
ing phase reveals the existence of two regimes. In the first 



w 



(5) 



where ^ is the correlation length, and Xm is a new expo- 
nent that characterizes the dependence of the prefactor 
of the width on the driving force. In the second regime, 
the effect of the pinned disorder becomes irrelevant com- 
pared to the annealed disorder [[16| , and we obtain 



r 



(6) 



where Xa is a new exponent, and aa is the roughness 
exponent corresponding to annealed disorder. Depending 
on the absence or presence of nonlinear terms, we recover 
the results of either the EW or the KPZ equation with 
annealed disorder. This crossover was observed both for 
experiments and simulations of discrete models . 

Since ^ ^ f~'^ , where v is the correlation length expo- 
nent, we propose the scaling ansatz 



(7) 



Upon comparison with (|^) and (^), we find that the 
scaling function g(u) satisfies g{u 3> 1) = const, and 
g{u ^ 1) m"'"^"''. We also obtain 



(8) 



In Fig. H we show the data collapse obtained using the 
results of Fig. ^^a) according to . The deviation from 
scaling for large values of £/^ is due to finite-size effects 
(see Fig. ||(b) for discussion). 

To determine a second relation for the new exponent 
Xm, let us consider the approach of the depinning transi- 
tion for a system of size L. For such a system the transi- 
tion does not occur for / = 0, as it would for an infinite 
system, but for an effective critical force such that f ~ 
implying / ~ L~^^''. Replacing this result into (^, we 
obtain 



£"" L°'^ 



(10) 



At the transition, (||) and ( |lO| ) must be identical. So, we 
find ap = a,„, and 



Xm = i' {ctL - dp). 



(11) 



Replacing the measured values of z/, a^, ap, and Ua into 
(|) and (|l]), we find Xm ^ 0.42 and Xa ^ 1-04, in good 
agreement with the values obtained in our simulations 
(Table |). 

Let us now consider the DPD universality class. For a 
description of the model, we refer to the literature [^,0. 
In Fig. 0(b), we show the local width for the pinned and 
the moving phases. We find ap ~ 0.63 for the pinned 
phase, and am — 0.75 for the moving phase (Table |). 
While previous studies considered am as an effective ex- 
ponent, and concluded that no self-afhne scaling exists 
in the moving phase for £ ^ ^ a consistent 

value for am regardless of how closely we approach the 
transition. This leads us to suggest the validity of the 
self-afHne scaling for the moving phase. This conclusion 
is supported by the good data collapse obtained with (|^), 
as shown in Fig. |^. 

An argument similar to the one leading to (O) yields 



V (ap 



(12) 
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Replacing the known values of v, ap, am, and aa into (|^) 
and (|l|), we find Xm — 

-0.2 and Xa ^ 0.23. Although 
the agreement with the measured values is not exact, the 
error bars do not rule out the validity of ([l^) . 

For the DPD universality class, we also find that /3 
changes values at the depinning transition (see Table |). 
We obtain for both sides of the transition a ~ /3, imply- 
ing that z remains unchanged at the transition jlj] . The 
exponent z characterizes the time scale tx for the propa- 
gation of correlations in the interface. This time scale is 
not expected to depend on the external force ||l^. These 
results are in good agreement with a numerical Integra- 
tionofEq. (0) @. 

Next, we compare our numerical results to earlier sim- 
ulations and discuss their significance for the determina- 
tion of the universality classes of several experiments. A 
problem with the interpretation of experimental and nu- 
merical results for the roughness exponent, has been the 
wide range of values for a: 0.5 — 1.0, measured in the mov- 
ing phase. We note that, in this regime, the crossover to 
the annealed disorder regime leads to effective exponents 
that change with the velocity (or the driving force). For 
this reason, we suggest that the scaling function (0) might 
be useful in the determination of the exponents from the 
study of the local width w. As shown in Table |, the 
exponent Xm has different signs for the two universality 
classes, leading to sharply distinct scaling behaviors for 
the prefactor of the width with the driving force. Since in 
many experiments it is possible to monitor the velocity of 
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the interface, and therefore the driving force, the study 
of this prefactor may also lead to an easier identification 
of the universality class to which the experimental results 
belong. 

In light of this discussion, the interpretation of the re- 
sults of Ref is clear. The numerical integration of the 
EW equation with quenched disorder performed in Ref. 
|pO| must belong to the QEW universality class, and the 
reported exponents are effective exponents whose values 
were affected by the crossover to the annealed regime. 

The determination of the universality class of the fluid 
invasion model (FIM) of Ref. is not trivial. However, 
we note that v is identical for the FIM and for the Hamil- 
tonian model. Furthermore, the value of a measured in 
Ref. might be explained as an artifact of the use of 
least square fits to the local width. 

Since the FIM is a realistic model for the fluid-fluid 
displacement experiments of Refs. we propose that 
they belong to the universality class of Eq. (||). In fact, 
the range of roughness exponents measured in the exper- 
iments of Refs. Q is consistent with the values that could 
be obtained with the Hamiltonian model. Also, in these 
experiments, it was found that the global width decreases 
with the increase of v, as for the QEW universality class. 

The paper burning experiments of Ref. |^ and the im- 
bibition experiments of Ref. Q are believed to belong to 
the DPD universality class. This conclusion is supported 
by the values of the exponents measured in both experi- 
ments. 

In summary, we find that for the DPD universality 
class, a and (3 change values at the depinning transition. 
We also find a new exponent Xm, that relates the values 
of the roughness exponents on both sides of the tran- 
sition. For the QEW universality class, a and /3 remain 
unchanged, and the exponent Xm relates the different val- 
ues of the local, ap, and global, ai,, roughness exponents 
at the depinning transition. 
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FIG. 1. (a) QEW universality class. Plot of the local width 
u) as a function of £ for several values of the driving force. The 
results were averaged over 30 realizations of the disorder and 
the system size is 5524. (b) DPD universality class. Plot of 
the local width ui as a function of £ for several values of the 
driving force. The system size is 10480 and each result was 
averaged over 50 realizations of the disorder. For clarity, the 
data for the pinned case is shifted by a factor of 1/2. 

FIG. 2. QEW umversaUty class. Data collapse, according 
to (7), of the results displayed in Figl(a). 
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FIG. 3. DPD universality class, (a) Data collapse accord- 
ing to (7), of the results displayed in Figl(b). (b) Data collapse 
for / = 0.0583 and systems of different size. It is visually ap- 
parent that the deviation from scaling observed in (a) is due 
to finite size effects. 

TABLE I. Critical exponents for the two universality 
classes studied in this paper. The subscript "p" refers to the 
pinned phase, the subscript "m" denotes the quenched disor- 
der regime in the moving phase, and the subscript "a" refers 
to the annealed disorder regime. The exponents z and Za 
were calculated from scaling relations, while the remaining 
exponents were calculated directly in the simulations. For the 
QEW, we obtain ul = 1.23 ± 0.04. 



Exponents 


DPD 


QEW 




0.63 ± 0.03 


0.92 ± 0.04 


Oim 


0.75 ±0.04 


0.92 ± 0.04 


Ota 


0.50 ± 0.04 


0.46 ± 0.04 


Z 


1.01 ±0.10 


1.45 ± 0.07 


Za 


1.67 ±0.26 


2.09 ± 0.40 


V 


1.73 ±0.04 


1.35 ± 0.04 


e 


0.64 ±0.12 


0.24 ± 0.03 


Pp 


0.67 ±0.05 


0.85 ± 0.03 


Pm 


0.74 ± 0.06 


0.86 ± 0.03 


Pa 


0.30 ±0.04 


0.22 ± 0.04 


Xm 


-0.12 ±0.06 


0.44 ± 0.05 


Xa 


0.34 ± 0.06 


0.99 ± 0.05 
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